arXiv:1507.04672v5 [hep-ph] 17 Jan 2017 


Axion-like particle dark matter in the linear regime of structure 

formation 


Qiaoli Yang^’^’* and Haoran Di^d 

^School of Physics, Huazhong University of Science and Technology, Wuhan 430074, China 
“^Department of Physics, Jinan University, Guangzhou 510632, China 

(Dated: January 18, 2017) 

Abstract 

If axion-like particles (ALPs) constitute a major part of dark matter (DM), due to the bosonic 
nature and a relative small mass, they could behave differently from the point like dark matter 
particles on the formation of the cosmic structures. When studying the structure formation, it is 
often useful to consider DM as a special fluid with a given density and a given velocity. ALP fluid 
obeys a same continuity equation comparing to the point-like collisionless DM but has a different 
first-order velocity equation. In the linear regime of structure formation, the resulted observational 
differences are negligible for the QCD axions but can be interesting for very light ALPs. 
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I. INTRODUCTION 


Q 


ain the strong CP problem 
3|. Generally speaking, a U(l) 


The introduction of the Peccei-Quinn mechanism to exp 

r 

leads to a new pseudoscalar particle known as the axion p, 
symmetry is introduced and spontaneously breaks down. A Goldstone boson appears and ac¬ 
quires a small mass from non-perturbative QCD effects thus it becomes a pseudo-Goldstone 
boson, the axion. Early axion models are quickly ruled out because their symmetry-breaking 
energy scales are not high enough to suppress their couplings to evade experimental detec¬ 
tions. Subsequent models such as the KSVZ axions and the DFSZ axions jd-T] have not 
been verihed due to the very weak interactions between the axions and the standard model 
oarticles. Thus, these models are called the ’’invisible axions.” Axion-like particles (ALPs) 
ilsi generically exist in the string theory, where antisymmetric tensor helds are compacti- 
hed on closed cycles, and the respective Kaluza-Klein zero modes acquire non-zero potentials 
due to non-perturbative effects on the cycles 10|[n|j^|. The potentials give a small mass 
to the Kaluza-Klein modes and these Kaluza-Klein modes are called the ALPs since they 
have many properties similar to that of the axions. The parameter spaces of these ALPs, 
however, are much less constrained than that of the QCD axions 13| 141. 

Cosmological and astronomical observations indicate that the majority of matter in our 
universe is composed of unknown particle types. From observations such as the Cosmic 
Microwave Background (CMB) anisotropy and primordial nuclear synthesis, we know that 
these particles are non-baryonic, cold, weakly interacting and stable. The other properties of 
these dark matter particles remain unknown. Axions/ALPs produced by the mis-alignment 
mechanism and/or decay of cosmic topological defects can have the required properties to 


represent a substantial fraction of dark matter 


3 - 201 . 


Currently, researchers are actively searching for cosmological or astrophysical original 
axions/ALPs, including ADMX 21| (Univ. of Washington), CAST 22| (Cern et al.). In 
contrast, PVLAS 23| (Padova) and ALPS 241 (DESY et al.) are aiming for the laboratory 
original axions/ALPs. 

Although the major candidates for dark matter, such as axions/ALPs, WIMPs and sterile 
neutrinos et al., were thought to be indistinguishable by purely astronomical and cosmologi¬ 
cal observations, the axions/ALPs have two interesting characteristics: they are bosonic and 
are highly occupied in the phase space. In addition, when treated as a special fluid, the hrst 
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order velocity equation of ALPs differs from that of the point-like dark matter particles, 
with additional terms that represent the quantum pressure and the self-interaction. The 
differences are negligible for the QCD axions in the linear regime, however, it is interesting 
to hnd whether these terms could have measurable consequences for extremely light ALPs 
with mass m 10-22eV. 


II. DARK MATTER AXIONS/ALPS 


The axion is a pseudo-Goldstone boson with a small mass due to the QCD instanton 
effect. The mass of an axion is determined by the potential term V{a) = A^[l — cos{a/fa)], 
where a is the axion held, fa is the symmetry-breaking scale and A is the QCD scale factor 
of order O.lCeV. From the axion potential, one can hnd the axion mass and self-interaction 
coupling by expanding near vacuum: m = A'^/fa, Cint = “(A"^a^)/(4!/^). Let us denote the 
interaction term as A = — A^/ ff for convenience. It is clear that the value of the symmetry¬ 
breaking scale fa is critical for axion phenomenology. The cosmological, astrophysical and 
experimental constraints on the value of fa are represented by lO^CeV < fa < lO^^CeV 


29 


35|. 


The axions are not absolutely stable because, for example, they can decay into two 
photons, but their decay rate is much smaller than the Hubble rate. Thus, practically 
speaking, the axions are stable. For a discussion of the properties of dark matter axions, 
one may use the following Lagrangian density: 


C 


^ 22 '^4 

-(oa) - m a —-a 

2^ ^ 2 4! 


( 1 ) 


The cosmology axions have two population types. The thermal axions are produced by 
thermal equilibrium with the primordial soup, so that the number density is determined by 
the decoupling temperature T^. Just before decoupling, we have [20|: 


n„ = 


1.202 3 Rz?, 
R ' 


TT^ 


( 2 ) 


where Rd and R are the scale factors for decoupling and for today, respectively. The 
decoupling temperature is determined by the condition T > H, where F is the rate axions 
are created and annihilated. The thermal axion density is very small compared to the critical 
density. 
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The second population type is the cold axions. The cold axions can be created by two 
ways: vacuum realignment and the decay of axion-held topological defects, such as the axion 
strings and/or the axion domain walls if the inflation happens before the PQ symmetry 
breaking. If the inflation happens after the PQ symmetry breaking, the topological defects 
are homogenized so that only the vacuum realignment contributes. For both scenarios, the 
cold axions have an abundance of order 20| |: 




fa 




( 3 ) 


TOi^GeV ’ 

which is a substantial fraction of the critical density if fa IQi^GeV. 

The ALPs are created due to compactifying antisymmetric tensor fields on closed cycles, 
and the respective Kaluza-Klein zero modes acquire a potential from non-perturbative effects 
on the cycles. There are many different ALPs due to the topology of the space-time manifold 
in the string theory and let us consider a scenario that one of them constitutes a substantial 
fraction of dark matter. The effective Lagrangian in the four dimensional space-time is: 

flip, 


C = 


-{da) - A^LpU{a) + .. , 


( 4 ) 


where /alp and Aalp are the decay constant and the potential energy scale for the ALPs, 
respectively. For small vibration of a we can always expand U around the minimum, so 
that the mass of the ALPs is of order m ~ A\^p/f alp- The potential energy scale Aalp 
depends on the UV energy scale of the string theory and also is exponentially sensitive to 
string instantons Sj. Therefore, the mass range of the ALPs is distributed on a large scale. 

The ALPs DM can be created in the early universe in a similar way as the QGD axions. 
The equation of motion of a homogenized ALPs held in a hat Friedmann-Robertson-Walker 
(FRW) universe is: 


d^a -|- ?>Ha -|- mfa = 0 , 


( 5 ) 


where H is the Hubble rate. After the Hubble rate decreases to comparing to the mass of 
the ALPs at time fi, ti ~ 1/m, the ALPs held starts to oscillate and the initial oscillation 
amplitude is determined by the vacuum realignment angle. The newly created ALPs are 
very cold so the energy density decreases as p oc R~^. The energy density of ALPs today 
is therefore p ~ {1 /2)m\fl[R{ti)/R{t)]^ and the ratio of the ALPs’ energy density to the 
critical energy density pc today is: p/p^ ~ < 1, where te is the matter-radiation 

equality time. In addition, the ALPs DM should be stable, which requires that the decay 
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rate: T ^ / (47r^/^) < Hq, where Hq is the Hubble rate today and a is the fine structure 

constant. The above two considerations give an upper limit of the mass of the DM ALPs: 
‘iixGvnJj < Hq, which implies m < 30keV. Finally, as the Compton wavelength 
of the DM ALPs should be smaller than Ikpc to allow the structure formation on the kpc 
scale, the DM ALPs have a lower limit on the mass m > 10“^^eV. 


III. AXION/ALP DM COMPARING TO THE POINT-LIKE COLLISIONLESS 
DM 


After the cold axions/ALPs are created by the vacuum realignment mechanism, these 
dark matter particles have a very high number density and a very small velocity dispersion. 
The high number density combined with a small velocity dispersion results in a very high 
phase space density. The phase space density Af of the QCD axion DM is |26| : 

fa , 


M ^ I0®i( 


10i2GeV^ 


( 6 ) 


For the ALPs lighter than the QCD axions, the phase space density is even higher. 

Because the quantum states is highly occupied by particles, one may use approxima¬ 
tion where operators become classical fields. Let us consider the axions/ALPs’ field in an 
unperturbed fiat FRW universe, the Lagrangian density is: 

C = - ^a^R^ + .. . (7) 


The equation of motion is thus: 


d‘fa — ^'V‘^a + 3Hdta + m‘^a + —a^ = 0 , ( 8 ) 

R^ 6 

where H = R/R is the Hubble parameter. 

The cold axions/ALPs are non-relativistic, so the slow varying term of the field is most 
important while the rapid varying terms of order axion/ALPs mass, such as can be 

factored out. Expressing that a = (-06“*”^* +'?/*e*™'*)/\/2m, where 0 is a complex scalar, the 
equation of motion for the slow varying term ijj is: 


— iip — 


1 


2mR? 


— i-Hijj + 


A 




V = 0 . 


(9) 
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To comparing the equation of motion with the point-like collisionless DM, it is useful to 
re-express Eq. ([9]) in terms of density and effective velocity. In comoving coordinates, one 
hnds the generalized continuity equation: 


dua , 1 d{naV^) , 

+ —— ?r- - h iHna = 0 , 


( 10 ) 


dt R dx^ 

where Ua = V'V'* is the particle density and u* = —i{'ip*di'ip — 'ipdi'ip*)/{2mRna) is the effective 
velocity. The continuity equation is the same as that of the point-like DM. However, the 
hrst-order velocity equation is different: 

dR „ , A „ 1 „ 


TIT + + 

at 




2m? R^ 


di 


n„ 


n„ 


= 0 , 


( 11 ) 


comparing with that of the point-like DM 


ot 


( 12 ) 


There are two additional terms: Xd^Ua/(Sm^R) and — (l/2m^i?^)c?j[(V The Erst 

one comes from the self-interaction and the second one comes from the quantum pressure. 

For the QCD axions, the self-interaction term and the quantum pressure term are too 
small to be considered for observable scales. This is obvious if one considers the Jeans length 
due to the self-interaction and the quantum pressure: 


di5 + 2Hdt5 + 


AttGp ■ 


r)^ = 0, 


(13) 


' Am^R^ R?' 

where 5 = 5p/p is the density contrast, p is the averaged density of DM and a is the cross 
section of the self-interaction. For QCD axions in the linear regime of density fluctuation, 
only very small-scale structures, {k/R)~^ rs-/ {y/a/Gm^Y^'^ and {k/R) ^ ~ {IQiim^Gp) 
differ from the point-like collisionless particle DM (see FIG.l). However, for the ALPs, the 
additional terms in Eq. flTT|) can be important if the mass is sufficiently small. 

In the non-linear regime of the evolution of axions/ALPs, the DM axions/ALPs thermal- 
ize due to the gravity interaction 2q, thus form a Bose-Einstein Condensation (BEC). If 
the wave function of the axions/ALPs can be written as: 


(14) 


and the effective velocity is defined as u = 9*0/ {Rm), the Landau’s postulates of a superffuid 
is then a natural consequence of the axions/ALPs BEC. In the non-linear regime of structure 
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FIG. 1: The quantum pressure versus the gravity (the quantum pressure dominates over the gravity 
below the red-line), and the self-interaction versus the gravity (the self-interaction dominates over 
the gravity below the black-dashed line) of the QCD axions. Thus in the linear regime of structure 
formation, the QCD axion DM behaves differently from the point-like collisionless particle DM 
only for very small scales ~ 10km. 


formation, superfluid behave very differently from point-like DM. For example, as V x u = 0, 
the angular momentum of the dark matter halos needs to be accommodated by appearance 
of quantized vortices. A detailed discussion of the formation of vortices in DM haloes is 
given in 


IV. SUMMARY 

We consider a scenario that the majority of dark matter is composed by axion like parti¬ 
cles. The cosmological and astrophysical considerations give rise to a mass constraint for the 
ALPs lO^^'^eV < m < 30keV. In addition, the hrst order velocity equation of the ALP fluid 
differs from that of the point-like collisionless DM particles with two additional terms that 
represent the quantum pressure and the self-interaction of ALPs respectively. The resulted 
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differences in the linear regime of structure formation are negligible for the QCD axions. 
However, it remains interesting to hnd whether these terms have measurable consequences 
for extremely light ALPs with mass m < 10“^^eV. 
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